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We show that, in contrast to classical random graph models, many real- world complex systems - 
including a variety of biological regulatory networks and technological networks such as the internet 
- spontaneously self-organize to a richly symmetric state. We consider the organizational origins of 
symmetry and find that growth with preferential attachment confers symmetry in highly branched 
networks. We deconstruct the automorphism group of some real-world networks and find that 
some, but not all, real-world symmetry can be accounted for by branching. We also uncover an 
intriguing correspondence between the size of the automorphism group of growing random trees and 
the random Fibonacci sequences. 

PACS numbers: 89.75.-k 89.75.Fb 05.40.-a 02.20.-a 



INTRODUCTION 

The use of coniplex networks to model the underly- 
ing topology of real- world complex systems has attracted 
much research interest [3, 0. Formally, the study of com- 
plex networks is underpinned by the mathematical the- 
ory of graphs and it is a classical result of random 
graph theory that almost all large random (Erdos-Renyi) 
graphs are asymmetric 3] . In this letter we consider the 
symmetry structure of a variety of real-world networks 
and find that, in contrast to the classical random graph 
models, a rich degree of symmetry is ubiquitous in real- 
world complex systems. 

Symmetry is a source of great abstract beauty'^ and is 
of great practical use in simplifying complex problems ^ . 
Accordingly, symmetry is a well-established cornerstone 
of many branches of physics 'g','^. However, despite a well 
developed abstract theory of graph symmetry]^, 0, ^3 ' 
the study of symmetry in real-world complex networks 
has been limited. 

In fact, when present, symmetry provides a pow- 
erful tool to understand network structure and func- 
tion. For example, symmetry can strongly affect a net- 
work's eigenvalue spectrum [ill - by giving rise to 
degenerate eigenvalues, for instance^] - causing it to 
deviate from those of the classical ensembles of ran- 
dom matrix theorv [l^ Il4j . Furthermore, symmetry can 
also have a profound influence on the dynamics of pro- 
cesses taking place on networks and may be ex- 
ploited to reduce the computational complexity of net- 
work algorithms [lil ll^ . 

Mathematically, a network is a graph, G = G{V,E), 
with vertex set, V (of size A^), and edge set, E (of size 
M) where 2 vertices are connected if there is an edge 
between them. Here we consider the symmetry struc- 
ture of various real-world networks via their automor- 



phism group. An automorphism is a permutation of the 
vertices of the network which preserves adjacency. The 
set of automorphisms under composition forms a group, 
Tq, of size agjl^j. If the network is a multi-digraph, 
we remove weights and directions and consider the au- 
tomorphism group of the underlying graph. A network 
is said to be symmetric (respectively asymmetric) if its 
underlying graph has a nontrivial (respectively trivial) 
automorphism group and the degree of network symme- 
try is quantified by ac- In order to compare networks of 
different sizes with each other we also consider the quan- 
tity tg — [ac/Niy^-^ which measures symmetry relative 
to maximum possible symmetry (the complete graph on 
N vertices and its complement, the empty graph on N 
vertices, are the most symmetric graphs, both having 
ac = N\). Other measures of graph symmetry are con- 
sidered elsewhere 0, lis]. Here, the nauty program 20] - 
which includes one of the most efficient graph isomor- 
phism algorithms |2l|| - is used to calculate the size and 
structure of the various automorphism groups. 

Table |I| gives the size of the automorphism group of 
some real- world complex networks, all of which are highly 
symmetric. Since classical Erdos-Renyi random graphs 
are generally asymmetric this rich degree of symmetry 
is surprising and begs an explanation. The ubiquity of 
symmetry in disparate real-world systems suggests that 
it may be related to generic self-organizational principles. 
In order to begin to investigate the relationship between 
symmetry and self-organization we consider how the pro- 
cesses of growth and preferential attachment affect sys- 
tem symmetry. 



SYMMETRY AND SELF-ORGANIZATION 

Based upon the observation that many real-world net- 
works are continuously growing and that new vertices 
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Network 


N 


M 


aa 


ra 


Human B Cell Genetic Interactions [22] 


5,930 


64, 645 


5.9374 X 10" 


4.6044 X 10"* 


C. elegans Genetic Interactions [23] 


2,060 


18,000 


6.9985 X 10^**^ 


1.5776 X 10"=* 


BioGRID datasetsf24l: 










Human 


7,013 


20, 587 


1.2607 X lO"**^ 


4.5418 X 10"" 


S. cerevisiae 


5,295 


50, 723 


6.8622 X lO*^* 


5.2753 X lO""* 


Drosophila 


7,371 


25, 043 


3.0687 X 10"^^ 


4.2993 X 10"" 


Mus musculus 


209 


393 


5.3481 X lO^^'^ 


5.1081 X 10"^ 


Internet at the Autonomous Systems Level[2^ 


22,332 


45, 392 


1.2822 X lo"'^''* 


3.9009 X 10"'' 


US Power Grid [26] 


4,941 


6,594 


5.1851 X lO^'^^ 


5.9011 X 10"'' 



TABLE I: The size of the automorphism group of some real-world networks. The size of the automorphism group of 
the giant component is given (to 5 significant figures) which, in all cases, contains at least 93% of the vertices in the network. 
Giant components were extracted using Pa1ek,27.l. 



often show a preference for attachment to more highly 
connected vertices, Barabasi and Albert proposed a sim- 
ple model which accounted for the origin of the power- 
law vertex degree distribution often seen in many real- 
world networks ;28j. Power-law degree distributions have 
been widely discussed 29] and a number of variations of 
the Barabasi- Albert model have been suggested 30, 
Here, we consider a simple variation based upon that 
proposed by Pennock et al^^. We start with the com- 
plete graph on k vertices as an initial seed. At each time 
step the system is updated by introducing a single new 
vertex and k new edges which connect the new vertex 
to those already in the system without allowing multiple 
edges. The probability that vertex Vi is chosen at time t 
to attach the new vertex to is given by 



„/ N dAt) , , 1 
P(w,;,t) = a ;/ + (1 - a) 



2M{t) 



Nit) 



(1) 



where M(i) = k{k— I) /2 + kt is the total number of edges 
in the system at time t, di{t) is the degree of vertex Vi at 
time t, Nit) = k + t\s the total number of vertices in the 
system at time t and < a < 1. Thus, when a = new 
vertices are attached to old in a purely uniformly random 
way, and when a — 1 new vertices are attached to old in 
a purely preferential way. 

Two distinct patterns of behavior emerge dependent 
upon k, the number of new edges associated with each 
new vertex. When k = 1, qq grows exponentially and 
the system becomes increasingly symmetric (see Fig^); 
while for k > 2, aa rapidly falls to zero, and the system 
becomes asymmetric for large time (see Fig^i). This 
discrepancy results since the two cases (fc = 1 and k > 2) 
represent two fundamentally different modes of growth. 
When k = I every new vertex is attached to the network 
with only one new edge so cycles cannot arise and the 
network is always a tree. Alternatively, when k > 2 new 
cycles are generated with each new vertex addition, and 
the network is never a tree. 

Parenthetically, these observations also suggest an in- 
triguing relationship between the size of the automor- 





FIG. 1: The effect of preferential attachment on net- 
work symmetry, (a) Growth of the automorphism group 
of growing random trees (k=l). Median trajectories over 50 
Monte-Carlo simulations for each a from — 1 in increments 
of 0.1 are shown (arrow shows increasing a). For all a the 
growth of the automorphism group shows a clear exponen- 
tial trend, (b) When fe = 1, limjv^cx) o-c'^ ~ '^i'^) ^'^^ ^ 
(median trajectories over 50 simulations are shown) . The top 
dotted line is y = (1.132-1-0.618)^ and the bottom dotted line 
is y = 1.132^. (c) When fe = 1, <j{a) grows exponentially with 
a showing that, for 'tree-like' growth, preferential attachment 
increases system symmetry, (d) Growth of the automorphism 
group for k > 2. Median trajectories over 50 simulations for 
fc = 5, 10, 15, 20, 25 are shown (arrow gives increasing k). Full 
lines give trajectories for a = 1, dotted lines give trajecto- 
ries for a = 0. When fc > 2 the networks quickly become 
asymmetric, with preferential attachment having little effect 
on system symmetry. 



phism group in growing random trees and the random 
Fibonacci sequences. The familiar Fibonacci sequence is 
generated by taking = 1, /2 = 1 and fi = -I- /^_2 
for i > 3. Similarly, the random Fibonacci sequences 
are generated by taking ri = 1, ^2 = 1 and = 
±ri_i ± ri-2 for all i > 3 where each ± is chosen inde- 
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FIG. 2: Growth of tree automorphism groups and the 
random Fibonacci sequences. In both panels, fc = 1. (a) 
The median trajectory for a — (full line) averaged over 
50 simulations and the line y = (1.132)^ (dotted), (b) The 
median trajectory for a = 1 (full line) averaged over 50 sim- 
ulations and the line y = (1.132 + 0.618)"= (dotted). 

pcndcntly with probability 0.5. Viswanath showed that 
hnii^ookiT/' = V = 1.13198824. . .:32j. Interestingly, 
our numerical simulations suggest that when k — 1, 
limAT^oo clq^ — <7{a) for all a where cr{a) is a con- 
stant (the Lyapunov constant js^) with cr(0) = v and 
a{l) = V + where = 1/2(1 -I- is the golden 
ratio (see Fig. ^ and Fig. 

Although the majority of large random Erdos-Renyi 
graphs are asymmetric |3, it is common for large ran- 
dom trees to exhibit a high degree of symmetry |3J|. In- 
tuitively, this is because the absence of cycles in trees 
means that the number of tree configurations are combi- 
natorially restricted (the set of trees on N vertices is a 
thin subset of the set of graphs on TV vertices) and this 
restriction can force repetition of identical branches from 
the same fork, endowing the tree with symmetries. For 
example, every tree contains at least two vertices of de- 
gree l[35|. These vertices are called leaves and are the 
end points of branches. Whenever a branch ends in 2 or 
more leaves an automorphism naturally arises by permu- 
tation of these leaves while holding all other vertices in 
the tree fixed. 

Consequently, we conclude that the rich degree of sym- 
metry seen in the k — 1 case is due to the fact that these 
networks grow as trees, rather than the mode of attach- 
ment of new vertices per se. However, we find that, given 
that the network is growing as a tree, preferential at- 
tachment increases network symmetry (see Fig^-c). A 
heuristic explanation for this is that preferential attach- 
ment introduces a bias toward multiple short branches 
and, probabilistically, short branches are more likely to 
be repeated about the same fork than longer branches. 
To see this consider the following two limiting cases: (1) 
in which new vertices are always attached to the most 
highly connected vertex. Starting from a single vertex 
this gives a t-stav, St, at time t where one vertex has de- 
gree t, all others have degree 1 and there are t branches 
of length 1. In this case, ag^ = (t — 1)! for t > 3 and 
the tree is almost maximally symmetric. (2) in which 
new vertices are always attached to the least highly con- 



nected vertex. Starting with a single vertex this gives a 
path (that is, a single long branch), Pt, at time t with 2 
vertices of degree 1 and t — 2 vertices of degree 2. In this 
case, = 2 for all t > 1 and the graph is minimally 
symmetric. 

These results suggest that branching may be a source 
of symmetry in real-world systems. In order to see if 
this is the case we considered the orbit structure of 
some real- world networks. For each vertex in a graph 
G, the orbit of v under the action of Tq is the set 
Tg{v) ^ {jv e G : -f e FgIH. When \TGiv)\ = 1 
the vertex is fixed (or trivial) otherwise it is non-trivial. 
Since orbits are disjoint^ the symmetry structure of a 
network can be investigated by visualizing its adjacency 
matrix and applying different colors to different orbits. 
As illustrative examples. Fig. Ogives the adjacency ma- 
trices of the human B cell and c. elegans genetic regu- 
latory networks. Here, in order further clarify network 
structure, the rows and columns of the adjacency ma- 
trices are also sorted in descending order by eigenvector 
centrality (eigenvector centrality is a measure of vertex 
importance in a network, and is the basis of numerous 
network algorithmsQ). Fig. shows that the majority 
of automorphisms of the B cell network are permutations 
of leaves indicating that, as in random trees, branching 
is the dominant source of symmetry in this network. By 
contrast. Fig. shows that the symmetry structure 
of the c. elegans network is much more intricate, with 
symmetry arising both from branching (as evidenced by 
movement of the leaves) and elsewhere. For example, the 
66 most central vertices in the c. elegans network prefer- 
entially associate with each other to form a highly cyclic, 
highly symmetric core. Of these 66 vertices, half asso- 
ciate exclusively with other vertices in this core. Since 
these 33 vertices can be permuted amongst themselves 
without affecting network structure, they form an orbit of 
length 33 (the largest orbit in the network) . The induced 
subgraph on these 33 vertices is complete, and therefore 
has an automorphism group of size 33! = 8.6833 x 10^^ 
(5 s.f.) which constitutes a significant proportion of the 
symmetry of the network as a whole. The presence of 
this highly cyclic, highly symmetric region demonstrates 
that not all real- world symmetry derives from branching. 

CONCLUSIONS 

In summary, we have shown that many real- world net- 
works are richly symmetric and have traced the origin of 
some of this symmetry to branching. Furthermore, using 
a simple mathematical model of network growth we have 
shown that preferential attachment can increase symme- 
try in tree-like regions by encouraging the formation of 
multiple short branches. However, as illustrated by the 
c. elegans genetic regulatory network, not all real-world 
network symmetry can be accounted for by branching 
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FIG. 3: Adjacency matrices of genetic interaction net- 
works sorted by eigenvector centrality. Edges associated with 
fixed leaves are darli gray crosses, tiiose associated witii non- 
trivial vertices are black circles and those associated with non- 
leaf fixed vertices are light gray points. Edges linking fixed to 
non-trivial vertices are colored as non-trivial, (a) The Human 
B cell genetic network, (b) The c. elegans genetic network. 



alone. 

Some potential benefits of grapli symmetry are known 
(symmetry enhances tolerance to attack 19], for exam- 
ple). However, in real- world networks it is unclear 
whether such benefits derive from underlying organiza- 
tional principles - with symmetry being incidental - or 
whether they are intrinsic to symmetry itself (growth 
with preferential attachment also enhances tolerance to 
attack 36], for example). In other words, it is unclear if 
symmetry is ever itself an organizing principle in net- 
work evolution. Elucidation of the relationship between 
symmetry and self-organization and further investigation 
of the symmetry structure of real-world networks is fer- 
tile ground for future research. We anticipate that such 
investigations will significantly enrich our understanding 
of the structure and behavior of many real- world complex 
systems. 
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